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ABSTRACT 


Recent  developments  and  applications  of  an  Interactive  boundary-layer  procedure 
for  unsteady  flows  are  reviewed.  The  emphasis  Is  on  a  model  problem  corres¬ 
ponding  to  an  oscillating  thin  airfoil  In  laminar  flows  and  results  are 
reported  for  different  amplitudes  and  frequencies  of  oscillation.  The  use  of 
the  characteristic  box  scheme,  with  Its  stability  criterion,  are  shown  to 
allow  the  accurate  calculation  of  reverse  flows  and  the  Interaction  procedure 
removes  the  singularity  to  allow  calculation  through  regions  of  separated  flow. 
Although  the  current  focus  of  the  Interactive  boundary-layer  procedure  has  been 
on  the  leading-edge  region,  It  has  general  applicability  and,  together  with 
models  for  transition  and  turbulent  flows.  It  can  provide  the  basis  for  a 
method  to  deal  with  oscillating  airfoils  and  wings  and  the  rapid  movement  of 
fixed-wing  arrangements  at  angles  of  attack  up  to  and  beyond  those  of  dynamic 
stall. 

Calculations  at  high  angles  of  attack  Indicate  that  the  oehavlor  of  the 
unsteady  separated  leading-edge  flow  has  similarities  to  steady  flows  down¬ 
stream  of  surface  corrugations .  The  use  of  linear  stability  theory  In  the 
latter  case  shows  that  the  location  of  the  onset  of  transition  moves  upstream 
with  severity  of  corrugation  and  can  move  Inside  the  separation  bubble.  In 
practice  this  means  that  the  bubbles  will  be  shortened  and  analogy  with 
unsteady  flows  suggests  that  transition  may  play  an  Important  role  and  pre¬ 
clude  the  existence  of  the  long  separation  bubbles  determined  by  the  laminar- 
flow  calculations. 
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1.0  INTRODUCTION 


The  experimental  Investigations  of  McCroskey,  McAlister  and  Carr  [1],  Carr, 
McAlister  and  McCroskey  [2],  Oadone  [3],  Francis,  Keese  and  Ret 1  lie  [4],  Daley 
and  Jumper  [5],  Francis  and  Keese  [6],  Lorber  and  Covert  [7]  and  Lorber  and 
Carta  [8]  were  prompted  by  the  need  to  provide  better  understanding  of  the 
flow  characteristics  of  oscillating  airfoils.  It  Is  Important,  for  example, 
that  the  rotor  blades  of  helicopters  do  not  experience  dynamic  stall  and  the 
experiments  have  demonstrated  that  stall  Is  preceded  by  a  near-surface  vortex 
associated  with  separation  of  the  unsteady  boundary  layer.  It  Is  likely  that 
the  vortex  stems  from  a  laminar  separation  bubble  which,  with  Increasing  angle 
of  attack,  grows  until  It  becomes  unstable  and  bursts  to  form  a  vortex  which 
rolls  and  grows.  These  experiments  have  been  complemented  by  a  number  of 
computational  Investigations  which  Involve  the  solution  of  the  Navler-Stokes 
equations  [9],  and  of  the  boundary-layer  equations  sometimes  coupled  to  Invls- 
cld  flow  equations  as,  for  example,  by  Cebecl  and  Carr  [10-12].  As  yet,  the 
calculation  methods  have  not  been  successful  In  representing  the  measurements 
and  a  procedure  to  provide  accurate  representation  of  oscillating  airfoil 
flows  for  a  wide  range  of  flow  conditions  Including  dynamic  stall  remains  to 
be  developed.  A  useful  review  covering  all  aspects  of  the  subject  has  been 
provided  by  Lang  [13]  and  Important  applications  considered  by  Herbst  [14]. 

The  calculation  of  flow  over  an  oscillating  airfoil  by  an  Interactive 
boundary-layer  scheme  makes  use  of  Invlscld  and  viscous  flow  methods  and 
couples  them  by  a  special  procedure.  The  success  and  accuracy  of  this  method 
depends  on  the  accuracy  of  each  component  Including  the  procedure  used  to  link 
the  Invlscld  and  viscous  flow  methods.  In  general  the  Invlscld  flow  methods 
are  based  on  panel  methods,  such  as  those  of  Gelssler  [15],  Maskew  and  Dvorak 
[16]  and  Teng  [17],  and  their  arcuracy  and  ability  to  solve  flow  problems  of 
great  complexity  are  well  known.  The  viscous-flow  methods,  on  the  other  hand, 
are  not  so  well  developed  and  their  application  to  time-dependent  flows  has  not 
received  the  same  Interest  as  to  steady  flows.  The  calculation  of  unsteady 
flows  with  large  regions  of  flow  reversal  has  bten  confined  mostly  to  model 
problems  Involving  flows  such  as  the  laminar  flow  over  a  circular  cylinder 
Impulsively  started  from  rest.  The  emphasis  has  been  on  the  solution  of  the 
boundary- layer  equations  for  a  prescribed  pressure  distribution  and  explora¬ 
tion  of  the  relationship  between  flow  singularity  and  separation.  In  spite  of 


the  apparent  simplicity  of  this  model  problem,  many  difficulties  were  encount¬ 
ered  and  some  remain  to  be  solved. 

The  calculation  of  boundary  layers  on  an  oscillating  airfoil  poses  some  prob¬ 
lems  that  are  different  than  those  considered  In  other  unsteady  boundary-layer 
flows.  Perhaps  the  most  Important  of  these  Is  the  calculation  of  the  upstream 
boundary  conditions  In  the  (t,y)  plane  at  some  x  *  xq  which  are  needed  In 
the  solution  procedure.  Contrary  to  steady  flows  where  the  streamwlse  veloc¬ 
ity  u  Is  Identically  zero  across  the  shear  layer  at  the  stagnation  point,  this 
Is  not  the  case  In  time-dependent  flows;  flow  reversals  occur  due  to  the  move¬ 
ment  of  the  stagnation  point  and  cause  the  locus  of  zero  u-veloclty  to  vary 
with  space  requiring  the  use  of  a  special  numerical  method.  It  should  be 
emphasized  that  the  problem  of  the  upstream  boundary  conditions  has  a  physical 
counterpart  and  Is  not  solely  numerical.  The  evolution  of  an  unsteady  flow  at 
high  angles  of  attack  can  Involve  the  formation  of  a  vortex  In  the  region 
close  to  the  leading  edge  and  Is  associated  with  the  onset  of  dynamic  stall. 
Thus  It  Is  Imperative  that  calculations  of  the  flow  In  the  upstream  region  be 
accurate. 

The  difficulties  In  the  solution  of  the  time-dependent  boundary-layer  equa¬ 
tions  depend  on  whether  or  not  there  Is  flow  reversal  across  the  boundary 
layer.  In  Its  absence,  there  are  several  numerical  methods  that  can  be  used 
to  solve  the  equations  Including  those  of  Crank-NIchol son  [18]  and  Keller  [19]. 
In  the  presence  of  flow  reversal,  however.  It  Is  necessary  to  use  a  scheme 
like  the  characteristic  box  method  which  Is  based  on  the  solution  of  time- 
dependent  boundary-layer  equations  along  the  local  streamlines  as  described  by 
Keller  [20]  and  Cebecl  [21]. 

This  scheme  allows  the  step  sizes  In  the  time  and  streamwlse  directions  to  be 
automatically  adjusted  to  ensure  that  the  region  of  backflow,  as  determined  by 
the  local  streamlines,  does  not  violate  a  stability  condition  like  that  of 
Courant,  Friedrichs  and  Lewy  (CFL).  Although  the  zig-zag  scheme  of  Krause  et 
al.  [22]  can  also  be  used  for  this  purpose,  It  can  be  Inaccurate  In  regions  of 
large  flow  reversal  since  the  orientation  of  numerical  mesh  Is  chosen  a  priori, 
as  discussed  by  Cebecl  [21]  and  Cebecl,  Khattab  and  Schlmke  [23]. 


The  coupling  between  Invlscld  and  viscous  flow  methods  has  attracted  even  less 
attention  than  the  above  and  most  work  has  again  been  directed  to  the  develop¬ 
ment  of  methods  for  steady  flows.  Two  coupling  procedures  have  been  proposed 
and  that  suggested  by  Veldman  [24]  offers  greater  flexibility  than  that  of 
LeBalleur  [25]  and  Carter  [26]  and  has  been  developed  for  two-dimensional  air¬ 
foils  by  Cebecl  and  his  associates  [27].  It  has  also  been  extended  to  three- 
dimensional  flows,  albeit  with  strip  theory  and  quasl-three-dlmenslonal 
approximations.  Including  flow  over  swept  wings  [28],  and  leading-edge  separa¬ 
tion  bubbles  on  thin  wings  [29]. 

The  flow  singularity  In  steady  flows  corresponds  to  the  vanishing  of  the  wall 
shear  and  the  boundary-layer  equations  can  be  solved  routinely  with  a  pre¬ 
scribed  pressure  distribution  up  to  this  point;  there  are  no  questions  regard¬ 
ing  the  accuracy  of  the  solutions.  The  u-veloclty  does  not  exhibit  flow 
reversals  across  the  layer  and  the  Inaccuracies  resulting  from  the  rapid  vari¬ 
ation  of  flow  conditions  near  the  separation  location  can  be  accounted  for  by 
taking  small  steps  In  the  streamwlse  direction.  It  Is  also  known  that  with 
Interaction,  that  Is  computing  the  external  velocity  as  part  of  the  solution, 
this  singularity  can  be  removed  and  separated  flows  can  be  calculated  for  both 
laminar  and  turbulent  flows.  Time-dependent  flows  follow  a  similar  pattern 
and  concentration  on  the  "proper"  definition  of  flow  singularity  or  separa¬ 
tion,  ensures  that  any  singularity  can  be  removed  by  Interaction.  This 
requires  that  the  singularity  Is  Identified  and  that  the  region  of  reverse 
flow,  which  usually  precedes  the  singularity  at  separation  of  an  unsteady 
flow,  be  calculated  accurately.  Inaccuracy  In  the  calculation  of  the  reverse 
flows  can  lead  to  breakdown  of  the  solution  procedure  and  the  erroneous 
assumption  that  It  Is  caused  by  a  singularity.  An  Interactive  boundary- layer 
method  may  not  remedy  this  situation  since  Its  success  depends  on  the  accuracy 
of  the  solutions  of  the  boundary-layer  equations.  We  shall  discuss  this  point 
further  In  this  report. 


2.0  MODEL  PROBLEM 


This  report  reviews  some  of  the  results  which  have  been  obtained  with  a  numer 
leal  procedure  designed  to  overcome  the  problems  outlined  above.  Thus,  with 
support  from  the  Air  Force  Office  of  Scientific  Research,  AFOSR ,  we  have 
addressed  fundamental  aspects  of  the  problem  Including  the  determination  of 
upstream  boundary  conditions,  the  relationship  between  singularity  and  sepa 
ration,  and  the  development  of  an  Interactive  procedure.  To  overcome  the 
corresponding  difficulties,  we  have  developed  a  novel  method  to  deal  with  the 
leading-edge  region  (Section  3.1),  an  accurate  numerical  method  for  the  region 
of  reverse  flow  (Section  3.2)  and  an  Interactive  procedure  which,  together 
with  a  numerical  stability  requirement,  have  allowed  us  to  study  the  nature  of 
the  singularity  and  to  provide  fundamental  understanding  of  the  mechanisms  of 
the  region  of  recirculating  flow  near  to  the  leading  edge  of  an  oscillating 
airfoil  (Section  3.3).  We  have  also  made  use  of  the  linear  stability  theory 
and  analogy  between  unsteady  flow  at  high  angle  of  attack  and  steady  flow  over 
a  surface  with  corrugation  to  conjecture  that  long  separation  bubbles  obtained 
for  laminar  flows  may  not  exist  In  practice.  It  appears  likely  that  transi¬ 
tion  will  occur  within  these  bubbles  with  consequent  reduction  In  their  length. 
This  Is  discussed  In  Section  4. 


To  accomplish  these  objectives,  we  have  considered  the  model  problem  of  lami¬ 
nar  flow  over  a  thin  airfoil  with  a  thickness  ratio  of  y  (=b/a)  at  a  reduced 

angle  of  attack  E  .  The  external  velocltv  for  a  steady  flow  In  this  case  can 
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be  t educed  from  Invlscld  flew  to  be 
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where  the  parameter  E  Is  related  to  the  surface  distance  s  by 
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Equation  (2) 
the  form 


s  =  ay2  /  (1  *  E2)1/2dE 
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Is  extended  to  unsteady  flows  by  Introducing  time  dependence  In 
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The  resulting  flow  contains  the  essential  Ingredients  of  the  leading  edge 
region  of  an  oscillating  airfoil  Including  the  moving  stagnation  point  witn 
consequent  reverse  flows,  the  possibility  of  boundary  layer  separation  and 
reattachment  and  their  dependence  upon  the  amplitude  of  oscillation,  A,  the 
frequency  w  and  distance  from  the  leading  edge,  s.  It  also  permits  compari 
son  with  the  equivalent  steady  flows  for  which  calculations  have  been  reported 
by,  for  example.  Cebecl,  Stewartson  and  Williams  [30]. 


The  coupling  between  the  Invlscld  and  viscous  flows  Is  accomplished  by  writing 

the  edge  boundary  condition  In  the  boundary- layer  equations  as  the  sum  of  the 
o 

Invlscld  velocity  t . t )  and  a  perturbation  velocity  Au  (E.t), 
that  Is, 
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and  Au^( £ , t )  Is  obtained  from  the  Hilbert  Integral  given  by 
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with  the  Interaction  region  confined  between  s^  and  s^  and  with  up(E.t)  given 
by  Eq.  (3).  Further  details  can  be  found  In  Ref.  30. 


The  emphasis  of  the  present  work  has  been  on  the  leading-edge  region  and  makes 
use  of  an  expression  for  the  freestream  boundary  condition  which  has  been 
obtained  by  solution  of  the  Invlscld-f low  equations.  The  Interactive 
boundary- layer  method  Is,  however,  general  so  that  It  can  Involve  any  solu 
tlons  to  Invlscld-f low  equations  and  an  entire  airfoil  with  laminar,  transl 
tlonal  and  turbulent  flow. 
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3.0  ACHIEVEMENTS 


The  following  three  subsections  highlight  results  obtained,  respectively,  with 
the  special  procedure  to  deal  with  the  stagnation  point,  to  show  the  Import¬ 
ance  of  numerical  accuracy  In  regions  of  reverse  flow  and  to  report  and  dis¬ 
cuss  the  effects  of  unsteadiness  on  recirculation  bubbles. 

3  1  Moving  Stagnation-Point  Problem 

The  numerical  difficulties  associated  with  the  moving  stagnation  point  exist 
in  a  number  of  flows  Including  those  associated  with  turbomachinery,  as  dis¬ 
cussed  by  Cebecl,  Slmoneau  and  Platzer  [31].  Thus  the  waxes  generated  Dy  mov¬ 
ing  rotor  blades  lead  to  movement  of  the  stagnation  points  on  the  downstream 
of  stator  blades.  An  accurate  numerical  scheme  has  been  developed  to  solve 
this  problem  and.  as  described  by  Cebecl  et  al.  [31],  makes  use  of  the  charac¬ 
teristic  box  scheme  in  an  iterative  manner  to  deal  with  the  flow  reversals. 


figures  1  and  2  show  calculated  results  for  the  model  problem  with  two  circu¬ 
lar  freguencles,  w  -  »/30,  »/3,  and  E  =  1.  A  =  -0.5.  They  allow  examination 

0 

of  the  effect  of  frequency  on  the  calculated  velocity  profiles  In  the  vicinity 
of  the  stagnation  point.  Figure  1  shows  that  the  locus  of  the  u-veloclty  on 
time  lines  t  =  t/2  and  *  Is  essentially  the  same  as  In  the  steady  case  and 
as  a  result  there  are  no  flow  reversals  In  the  velocity  profiles.  However,  as 
can  De  seen  from  Figure  2,  flow  reversals  begin  to  occur  around  the  stagnation 
point  with  increase  in  the  frequency  to  w/3  and  become  more  prolonged  as  time 
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Fig.  2.  Variation  of  the  velocity  profiles  for  w  =  »/3  near  the  stagna¬ 
tion  region  for  different  values  of  wt,  (a)  90,  (b)  180,  (c)  270, 

(d)  360.  The  dashed  line  Indicates  the  locus  of  zero  u-veloclty. 

Increases  to  wt  =  3*/2.  At  wt  =  2»,  the  region  of  flow  reversal  Is  reduced 
but  Is  not  zero  as  It  was  at  wt  =  0. 


3.2  Large  Flow  Reversals  and  the  Need  for  Numerical  Accurac 


Our  calculations  for  the  model  problem  were  arranged  to  parallel  those  reported 
by  Cebecl  [21]  for  the  flow  over  a  circular  cylinder.  Thus  the  zig-zag  and 
the  characteristic  box  schemes  were  used  with  time  and  distance  steps  chosen 
arbitrarily  and  the  calculations  repeated  with  values  In  agreement  with  a  sta¬ 
bility  criterion.  The  results  of  Fig.  3  for  w  =  0.10  were  obtained  with  the 
zig-zag  box  scheme  by  Cebecl,  Khattab  and  Schlmke  [32]  for  a  A£-spaclng 
specified  such  that  A^  =  0.01  up  to  5  =  1.7,  A^  =  0.005  for  1.7  <  £  <  4  and 

AE  =  0.01  for  4  <  l  <  8;  the  time  steps  At  were  10  degrees  for  0  <  wt  <  260°, 
1  n 

5  degrees  for  260°  <  wt  <  295°,  and  1.25  degrees  for  295°  <  wt  <  360°.  The 
calculations  broke  down  at  wt  =  310°,  suggesting  flow  separation  at  this 
location. 


Fig.  3  shows  that  the  variation  of  the  dimensionless  displacement  thickness 
6*  Is  generally  smooth  except  In  the  neighborhood  of  t  =  2.12  and  for  wt 
*  308.75°.  Here  R  Is  Reynolds  number  based  on  nose  radius  and  Is  defined  by 
(2au«/u).  The  first  sign  of  Irregularity  Is  the  steepening  of  the  slope  of 
of  4*  when  wt  =>  300°  and  a  local  maximum  of  4*  occurs  at  l  =  2.12  when  wt  = 
308.75°.  When  the  same  results  were  plotted  for  a  displacement  velocity,  (d/d£) 
(ue4°),  we  observed  steepening  of  the  displacement  velocity  as  the  peak 
move  from  £  =  2.125  to  2.08  with  wt  changing  from  300  to  308.75  degrees. 

The  calculations  which  led  to  Fig.  3  were  repeated  with  the  characteristic  box 
scheme  using  the  same  coarse  variations  of  At  and  A£.  and  the  results 
were  Identical  to  those  obtained  with  the  zig-zag  scheme  up  to  wt  =  280°.  At 
wt  =  282.5°,  the  solutions  of  the  zig-zag  scheme  were  smooth  and  free  of 
wiggles  but  those  of  the  characteristic  box  scheme  exhibited  oscillations  In 
f  which  led  to  their  breakdown.  The  solutions  with  the  zig-zag  scheme, 
however,  continued  without  numerical  difficulties  until  wt  =  310°,  where 
oscillations  appeared  and  led  to  the  breakdown  of  the  solutions  at  the  next 
time  step. 

The  character  1 stlc  box  was  used  subsequently  with  values  of  A^  fixed  as 

before  and  with  values  of  At  determined  In  accord  with  the  stability 

n 

requirement  as  discussed  by  Cebecl  et  al.  [23].  This  procedure  avoided  the 


Fig.  3.  Variation  of:  displacement  thickness,  4*  with  F  for  the  oscillat¬ 
ing  airfoil,  A  =  1,  w  =  0.1,  obtained  with  the  zig-zag  scheme. 
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breakdown  of  the  solutions  and,  as  can  be  seen  from  Fig.  4,  the  maximum  value 
of  the  stability  parameter  3  Increases  considerably  with  wt  so  that  the 
solutions  required  correspondingly  smaller  values  of  the  time  step.  It  Is 
Interesting  to  note  that  the  wall-shear  distributions  of  Fig.  5  are  uninflu¬ 
enced  by  the  mesh  at  wt  =  280#  and  310°  but,  for  wt  >  310°,  the  coarse 
mesh  leads  to  large  values  of  8  and  breakdown  of  the  solutions. 

Figure  6a  shows  the  distributions  of  displacement  thickness  for  values  of  wt 
from  260°  up  to  360°  and  completes  the  cycle.  The  results  up  to  300°  were 
Identical  with  those  of  Fig.  3  with  rapid  Increase  of  the  displacement  thick¬ 
ness  corresponding  to  the  Increasing  extent  of  flow  reversal,  as  shown  by  the 
wall-shear  distributions  of  Fig.  6b. 

It  can  also  be  seen  from  this  figure  that  the  maximum  displacement  thickness 
and  minimum  wall  shear  move  upstream  with  Increasing  wt  for  values  of  wt 
up  to  324.5°;  this  feature  was  also  observed  In  the  calculations  performed  for 
the  circular  cylinder  [21].  The  results  obtained  with  the  zig-zag  scheme  and 
values  of  flt^  determined  by  the  characteristic  scheme  for  the  oscillating 
airfoil  were  Identical  to  those  discussed  above,  and  similar  correspondence 
was  obtained  with  the  calculations  performed  for  the  circular  cylinder. 


Fig.  4.  Effect  of  the  coarse  and  fine  meshes  on  the  variation  of  the 
stability  parameter  8  with  wt . 
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Effect  of  the  coarse  and  fine  meshes  on  the  variation  of  the  (a) 
stability  parameter  0,  and  (b)  wall  shear  f£  with  E. 


6.  Results  obtained  with  the  characterlstlc^box  scheme  for  u  =  0.1. 
V^j1a£1on  of  (a)  displacement  thickness  i*  and  (b)  wall  shear 


3.3  Effects  of  Unsteadiness  and  Interaction  on  Separated  Flows 
Figure  7  shows  distributions  of  wall  shear  for  frequencies  w  =  0.01  and 
0.001.  As  expected,  the  critical  value  of  the  reduced  angle  which  corresponds 
to  separation.  Is  smaller  than  that  for  the  higher  frequency  of  Section  3.2 
and  closer  to  that  of  the  steady  state,  E  =  1.16.  For  u  =  0.01,  the 
breakdown  of  the  solutions  occurs  at  wt  =  226®,  which  corresponds  to  an 
effective  reduced  angle  of  Ee^  =  1.360;  for  «*>  =  0.001,  the  correspond¬ 
ing  values  are  wt  =  204®  and  1.203.  We  also  note  from  Figs.  7a  and  7b  that 
the  flow  Is  a  little  unsteady  even  at  these  frequencies,  and  the  solutions  do 
not  break  down  with  the  appearance  of  flow  reversal,  which  Increases  In  extent 
as  u  changes  from  0.001  and  0.01. 

The  above  results,  and  those  of  the  previous  subsection,  were  obtained  without 
Interaction  between  the  solutions  of  Invlscld-f low  and  boundary-layer  equa¬ 
tions.  In  contrast,  the  results  of  Fig.  8  were  obtained  with  Interaction  and 
were  performed  In  the  following  way.  For  all  values  of  time  with  wt  ranging 
from  0®  to  360®,  the  standard  method  and  the  leading-edge  region  procedure 
discussed  In  Section  3.1  were  used  to  generate  upstream  boundary  conditions  at 
a  short  distance  from  the  leading  edge,  E  =  0.5.  With  these  conditions  and 
for  each  value  of  wt,  the  Inverse  method  was  used  to  calculate  the  unsteady 
flow  from  E  =  0.5  to  10,  for  the  specified  value  of  Reynolds  number  and 
thickness  ratio  of  y  =  0.1.  Since  the  system  of  equations  Is  now  elliptic, 


(a)  (b) 


Fig.  7.  Variation  of  wall  shear  fw  with  E  for  (a)  u  =  0.01,  (b)  w  =  0.001. 


Fig.  8.  Effect  of  Interaction  on  the  variation  of  wall  shear  fw  with  E 
for  (a)  w  =  0.01,  (b)  u.  =  0.001.  R  «  10*. 

sweeps  In  the  E-dlrectlon  were  necessary  to  achieve  a  converged  solution; 
around  three  sweeps  were  required  where  flow  reversal  was  encountered  and  a 
single  sweep  sufficed  where  It  did  not.  It  Is  to  be  expected  that  the  value 
of  R  will  Influence  the  number  of  sweeps  and,  since  It  Is  linked  to  physical 
parameters,  will  affect  the  singularity  and  the  size  of  the  bubble. 

Figure  8  shows  the  results  for  u  =  0.01  and  0.001  with  R  =  10*.  They  are 

nearly  the  same  as  those  obtained  by  the  standard  method  as  shown  in  Fig.  7 

prior  to  flow  reversal  where  the  Influence  of  the  Reynolds  number  Is  small  and 

Increase  after  flow  reversal.  In  the  case  of  u  =  0.001,  for  example,  the 

standard  method  predicts  flow  reversal  around  E  =  1.19  (see  Fig.  7)  and 

ef  f 

with  Interaction  (Fig.  8)  this  effective  angle  Is  between  1.219  and  l  2b4 

M 

The  maximum  negative  value  of  the  wall  shear  parameter  f^  obtained  with  the 

standard  method  Is  around  -0  03  at  E  ,,  =  1.199  and  may  be  compared  with  the 
„  ef  f 

maximum  value  of  fy  of  -0  14  at  Ee^  =  1.286  obtained  with  interaction  As 
expected,  the  Interaction  allows  the  calculations  to  be  performed  at  higher 

angles  of  attack  than  those  achieved  with  the  standard  method.  For  w  .  0.001, 
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the  maximum  for  which  calculations  can  be  performed  with  the  standard 

method  Is  1.199  with  breakdown  occurring  at  £  =  1.209;  the  corresponding 

ef  f 

values  with  Interaction  are  1.286  and  1.287.  Comparison  of  wall-shear  results 
with  both  procedures  and  w  =  0.001  Indicates  that  the  extent  of  the  recircu¬ 
lation  region  A£  Is  around  0.5  for  the  standard  case,  and  around  2.5  for  the 
Interactive  case.  The  solutions  do  not  have  a  singularity  In  the  former  case 
but  do  contain  flow  reversals  and  this  suggests  that  time-dependent  flows  can 
be  calculated  without  using  an  Inverse  procedure.  As  the  angle  of  attack 
exceeds  £  ^  =  1.199  for  u  =  0.001,  a  singularity  develops  and  requires  an 
Inverse  procedure  as  In  two-dimensional  steady  flows.  This  procedure  allows 
the  calculation  of  larger  regions  of  reverse  flow  where  the  flow  Is  now 
separated. 

We  see  a  similar  picture  with  the  greater  unsteadiness  corresponding  to  w  = 

0.01,  for  which  the  standard  method  allows  calculations  up  to  an  effective 

angle  of  attack  of  1.354  (Fig.  7a),  a  value  considerably  higher  than  1.199 

obtained  at  w  =  0.001.  The  first  flow  reversal  occurs  shortly  after 

l  ,c  =  1.294  and  breakdown  occurs  at  £  ,,  =  1.360  with  maximum  negative  wall 
eff  eff 


shear  values  of  -0.14  at  £  ,,  =  1.354  and  -0.035  at  £ 

eff  eff 


1.315.  The  extent 


of  the  maximum  reverse-flow  region  Is  now  1.5,  considerably  larger  than  for 

w  =  0.001,  and  Indicates  that  the  more  unsteady  nature  of  the  flow  produces 

a  bigger  region  of  reverse  flow  free  from  singularities.  For  this  value  of 

w,  the  Interactive  scheme  Increases  the  value  of  £  ,,  for  which  solu- 

tlons  can  be  obtained  to  1.424  with  breakdown  occurring  shortly  after  this 

value  at  1.428  (see  Fig.  8b).  The  first  flow  reversal  occurs  after  £  rr  = 

eff 

1.315  with  maximum  negative  wall  shear  equal  to  -0.19  at  £  =  1.424, 

eff 

and  the  extent  of  the  recirculation  region  has  now  Increased  by  about  30%. 

II 

Comparison  of  maximum  wall  shear  values,  f  ,  at  a  similar  value  of  £  ,,  1 nd 1 

w  eff 

cates  that  those  computed  with  the  Interactive  scheme  are  lower  than  those 
with  the  standard  scheme  so  that,  for  example,  the  Interactive  scheme  gives 

M 

(f  )  =  -  0.04  at  £  ff  =  1.36  compared  to  -0.14  at  £  ff  =  1.354  with  the 

standard  method  (Fig.  7a). 


Since  the  calculations  began  at  wt  =  0  with  solutions  obtained  by  solving 
steady-state  equations,  It  was  necessary  to  confirm  the  extent  of  their  Influ 
ence.  As  a  consequence,  calculations  were  performed  for  a  second  cycle  and 
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gave  results  at  wt  »  720°  which  were  Identical  to  those  at  wt  =  360°,  con¬ 
firming  that  the  flow  Is  cyclic.  Examination  of  the  results  showed  that  the 
Influence  of  the  Initial  conditions  die  out  rapidly  and  have  no  Influence  on 
the  solutions  presented  here. 

The  results  obtained  with  u>  »  0.001  can  usefully  be  compared  with  the  steady- 
state  results  of  Cebecl  et  al.  [30]  shown  In  Fig.  9.  We  might  expect  that  the 
small  unsteadiness  associated  with  this  frequency  will  lead  to  results  very 
similar  to  those  of  steady  state.  Inspection  of  Figs.  7b  and  9  shows  that 
although  this  Is  correct  In  general  terms,  the  answers  are  quantitatively 
different.  As  can  be  seei  ,  the  maximum  effective  angle  at  which  solutions  can 
be  obtained  Is  greater  In  the  unsteady  case  by  some  7%.  There  are  differences 
In  the  two  calculation  procedures  but  It  Is  unlikely  that  they  are  responsible 
for  this  difference.  On  the  other  hand.  It  Is  possible  that  the  difference  In 
the  negative  wall  shear  values  may  have  been  Influenced  by  the  use  of  the  FLARE 
approximation  In  the  steady-state  solutions.  Nevertheless,  the  unsteady  nature 
of  the  flow  with  u  »  0.001  Is  clear.  In  spite  of  this  very  low  reduced  fre¬ 
quency. 


Fig.  9.  Effect  of  Interaction  on  the  variation  of  wall  shear  parameter 
f^  for  a  steady  flow  at  R  =  104. 


4.0  CONJECTURES;  SEPARATION  BUBBLES  ANO  TRANSITION 


An  Important  conclusion  from  the  results  of  Section  3.3  Is  that  the  problem  of 
the  singularity,  discussed  In  several  papers  Including  that  of  Cebecl  et  al. 
[32],  has  been  overcome  by  the  combined  use  of  an  Interactive  boundary-layer 
method  and  a  novel  numerical  scheme  with  Its  stability  requirements  fulfilled. 
This  combination  had  allowed  calculations  to  be  performed  but  there  remains  a 
concern  about  the  physical  realism  of  the  long  regions  of  reverse  flow  pre¬ 
dicted  In  some  circumstances.  Experiments  suggest  that  leading-edge  regions 
of  reverse  flow  tend  to  be  small  with  transition  to  turbulent  flow  making  dif¬ 
ficult  the  existence  of  long  bubbles.  This  concern  Is  examined  by  Cebecl  [33] 
In  the  context  of  steady  flows.  It  Is  briefly  described  In  the  following 
paragraphs  and  the  Implications  of  the  results  for  unsteady-flow  calculations 
considered. 

As  a  steady-flow  counterpart  of  the  above  flows,  calculations  were  performed 
for  the  thin  airfoil  as  a  function  of  angle  of  attack  and  for  the  same  thick¬ 
ness  ratio  y  (=0.1)  and  were  reported  by  Cebecl  et  al.  [30]  and  Cebecl  [33] 

6  5 

for  Reynolds  numbers  of  10  and  10  ,  respectively.  The  corresponding  dls- 

II 

trlbutlons  of  f  are  shown  on  Fig.  10  and  have  the  same  form  as  those 

w 

presented  for  the  oscillating  airfoil.  At  the  highest  angle,  the  separation 
bubble  tends  to  grow  with  the  number  of  numerical  Iterations  and  suggests  that 
the  bubble  Is  unstable. 

The  large  separation  bubbles  calculated  for  steady  flow  seem  unrealistic  and, 
at  least  for  the  10^  Reynolds  number  where  transition  Is  certain  to  occur 
downstream  of  the  bubble,  there  Is  the  possibility  that  the  real  flow  cannot 
remain  laminar  In  the  region  of  Interest  and  that  this  may  affect  the  bubble 
length.  There  are,  however,  no  experimental  data  for  airfoils  with  which  to 
test  the  possibility  and  recourse  has  been  made  to  a  different  flow  for  which 
Information  of  transition  Is  available  and  which  has  close  similarities  to  the 
leading-edge  flow. 

The  experiments  of  Fage  [34]  were  performed  with  a  series  of  bumps  on  a  plate 
and  provide  Information  of  the  location  of  transition  downstream  of  a  region 
of  separated  flow  caused  by  the  bump  and  the  consequent  adverse  pressure 


Fig.  10.  Separation  bubbles  and  their  breakdown  -  explanation  from  steady 
flow  (l.e.  w  =  0),  R  =  105. 

II 

gradient.  Figure  11  shows  the  bump  and  a  series  of  distributions  of  f 
calculated  with  the  Interactive  procedure:  more  extensive  results  have  been 

n 

reported  by  Cebecl  and  Egan  [35].  The  shapes  of  the  f  distributions  are 

w 

similar  In  form  to  those  of  Figs.  8  and  9,  downstream  from  the  beginning  of 
the  favorable  pressure  gradient  and  lend  support  to  the  Idea  that  conclusions 
based  on  the  bump  flows  can  be  extrapolated  to  those  on  steady  and  unsteady 
thin  al rf ol 1 s . 

The  calculations  which  led  to  the  results  of  Fig.  17  provided  velocity  pro¬ 
files  and  those,  In  turn,  were  used  In  the  en-method  of  Smith  and  Gamberonl 
[36]  and  of  Van  Ingen  [37]  to  calculate  the  location  of  transition  shown  on 
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Fig.  11.  Accuracy  of  stability  approach  to  transition.  Steady  flow  over 
bumps,  data  of  Fage  [34]. 

the  figure.  This  method  solves  the  Orr-Sommerfeld  equation  for  given  velocity 
profiles  and  computes  the  amplification  rates  to  determine  the  location  of  the 
onset  of  transition.  As  can  be  seen  from  Fig.  11,  the  calculated  values  of 
the  location  of  transition  agree  well  with  measured  values  and  Indicate  a  clear 
trend  for  transition  to  move  upstream  with  Increase  In  the  height  of  the  bump 
and,  therefore,  with  the  strength  of  the  adverse  pressure  gradient  and  length 
of  the  region  of  separation.  In  the  case  of  the  largest  bump  shown,  the  tran¬ 
sition  location  Is  Inside  the  bubble  and.  If  the  Interactive  boundary- layer 
calculations  were  to  be  repeated  with  transition  specified  at  this  location, 
the  bubble  would  be  correspondingly  reduced  In  length. 
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It  can  be  conjectured  that  the  separation  bubble  would  be  further  shortened  as 
the  adverse  pressure  gradient  Is  Increased  In  strength  and  transition  moves 
further  upstream.  This  result  can  be  extrapolated  to  the  steady  thln-alrfoll 
flows  discussed  above  and  their  unsteady  counterpart.  Thus,  It  seems  likely 
that  the  long  separation  bubbles  calculated  above  are  unrealistic  and  In 
practice  would  be  shortened  by  transition  from  laminar  to  turbulent  flow. 
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5.0  CONCLUDING  REMARKS 


A  series  of  contributions  have  been  made  to  the  understanding  of  the  flow 
characteristics  of  the  leading-edge  region  of  an  oscillating  airfoil.  They 
required  the  development  and  use  of  a  calculation  method  which  Involves  the 
accurate  solution  of  time-dependent  boundary-layer  equations  and  their  inter¬ 
action  with  solutions  to  Invlscld-f low  equations.  The  need  for  accurate  calc¬ 
ulations  of  regions  of  reverse  flow,  Including  those  associated  with  the  mov 
Ing  stagnation  point,  demanded  the  use  of  the  characteristic  box  scheme  and 
the  fulfillment  of  an  associated  stability  requirement.  Although  those  formu¬ 
lated  Investigations  have  been  made  In  relation  to  the  model  problem  of  laminar 
flow  In  the  leading-edge  region  of  a  thin  airfoil,  the  calculation  method  Is 
general  and  can  readily  be  extended  to  permit  the  calculation  of  the  flow 
around  an  oscillating  airfoil. 

Among  the  specific  conclusions  which  can  be  drawn  are  the  following: 

1.  The  accuracy  of  the  results  obtained  from  the  solution  of  the  boundary 
layer  equations  has  been  examined  with  emphasis  on  regions  of  flow  rever 
sal  and  separation  where  the  characteristic  box  scheme  Is  used.  Attempts 
to  Improve  accuracy  by  ad  hoc  changes  to  the  finite-difference  mesh  failed 
and  revealed  the  need  for  a  procedure  which  would  automatical ly  guarantee 
accuracy  by  the  selection  of  an  appropriate  mesh.  This  was  achieved 
through  a  stability  criterion,  similar  to  that  of  Courant,  Friedrichs  and 
Lewy  [38].  The  combination  of  this  requirement  and  the  characteristic 

box  scheme  led  to  accurate  solutions  and  showed  that  the  mesh  require¬ 
ments  were  extremely  severe  In  the  region  of  large  flow  reversals. 

2.  Calculations  performed  for  a  range  of  reduced  frequencies  from  0  to  0.1 
show  that  Increased  unsteadiness  allows  results  to  be  obtained  at  higher 
angles  of  attack  before  the  solutions  break  down  and  that.  In  the  case  of 
the  highest  frequency,  there  was  no  breakdown.  The  calculations  with  the 
standard  method  led  to  regions  of  flow  reversal  whlcn  were  limited  In 
their  extent  by  the  singularity  except  at  the  highest  frequency.  The 
Interactive  procedure  removed  this  singularity  and  resulted  In  larger 
regions  of  flow  reversal  which  Involved  separation  at  higher  angles  of 
attack . 


3. 


The  unsteady  nature  of  the  flow  at  the  highest  frequency  allowed  the 
calculation  of  large  regions  of  flow  reversal  and  It  Is  expected  that  yet 
higher  frequencies  will  lead  to  even  larger  regions  of  flow  reversal. 

This  In  turn  will  permit  calculations  to  be  performed  at  larger  angles  of 
attack  where  the  occurrence  of  the  singularity  will  require  the  use  of 
the  Interactive  procedure.  The  gains  In  angles  of  attack  are  again  likely 
to  be  limited  by  the  ability  of  the  laminar  flow  to  sustain  separation 
bubbles  and  results  obtained  from  the  application  of  linear  stability 
theory  suggest  that  the  location  of  the  onset  of  transition  moves  upstream 
with  Increasing  angle  of  attack  and  extent  of  the  separation  bubble  so 
that  It  can  occur  within  the  bubble.  In  practice,  this  would  lead  to  the 
curtailment  of  the  region  of  separated  flow  and  the  development  of  a  turb¬ 
ulent  boundary  layer.  Support  for  the  accuracy  of  the  calculations  of  the 
onset  of  transition  Is  afforded  by  comparison  of  calculated  results  with 
measurements  in  the  steady  flow  over  bumps  which  give  rise  to  distribu¬ 
tions  of  the  wall-shear  parameter  similar  to  those  of  airfoil  flows. 

The  research  which  led  to  the  preceding  conclusion  Is  fundamental  In  nature  but 
has  practical  implications  which  should  be  emphasized.  The  calculation  method 
has  been  applied  to  a  model  problem  and  shown  to  be  able  to  represent  unsteady 
laminar  flows  with  regions  of  flow  reversal  and  separation.  It  can  equally  be 
applied  to  calculate  the  flow  over  complete  airfoils  subjected  to  cyclic  or 
other  time-dependent  boundary  conditions.  The  additional  features  to  allow 
these  extensions  already  exist  so  that,  for  example,  turbulence  models  have 
been  embodied  In  a  steady- state  version  of  the  method  and  have  been  shown  to 
represent  boundary-layer  and  wake  flows  with  accuracy  close  to  that  of  the  test 
measurements;  steady-state  calculations  of  transitional  flows,  Including  those 
with  flow  separation,  have  been  examined  and  shown  to  be  calculable;  and  a 
foundation  has  been  laid  for  the  calculation  of  flows  around  components  of 
aircraft.  With  this  catalogue  of  achievements.  It  Is  recommended  that  work  be 
carried  out  to  develop  the  method  further  and  to  apply  It  to  problems  Includ¬ 
ing  those  associated  with  helicopter  blades  and  the  maneuverability  of  air¬ 
craft  such  as  the  X - 3 1  so  that  the  merits  of  different  designs  can  be  assessed. 
It  Is  particularly  appropriate  that  It  be  used  to  study  flows  associated  with 
angles  of  attack  close  to  those  of  static  and  dynamic  stall. 
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